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On  The  Classification  and  Nonlinear  Solution  of  One  Type  of  / 121 

Stratified  Shear  Flow 

Liu  Shida  Chen  Jiayi  Liu  Shishi  Wang  Shufang 
(Department  of  Geophysics,  Beijing  University) 


ABSTRACT 

Starting  from  the  nonlinear  internal  gravity  wave  equation 
for  stratified  shear  fluid,  this  paper  considered  a  type  of 
progressive  wave  motions  and  derived  a  set rof  two-variable 
(perturbation  velocity  and  perturbation  density)  ordinary 
differential  equations  for  the  self-governing  dynamic  system  of 
this  type  of  fluid  motion.  Based  on  qualitative  theories  in 
differential  equation,  this  paper  also  conducted  a  qualitative 
analysis  on  the  geometric  topological  structure  of  the 
integration  curve  in  the  vicinity  of  the  origin  of  a  phase  plane 
which  used  the  perturbation  velocity  and  perturbation  density  as 
coordinates.  Based  on  the  difference  in  the  Richardson  number  on 
the  phase  diagram,  the  integration  curve  was  separated  into 
several  areas  with  different  properties.  When  Ri<0,  the  singular 
point  was  an  unstable  saddle  point  regardless  of  whether  the 
velocity  shear  du/dz  was  positive  or  negative.  When  Ri>0,  the 
singular  point  was  unstable  in  an  area  where  the  velocity  shear 
du/dz  was  positive,  while  the  singular  point  was  stable  where  the 
velocity  shear  du/dz  was  negative  (when  0<Ri<1/4,  it  was  a  stable 
nodal  point;  when  Ri>1/4  it  was  a  stable  focal  point,  so  the 
periodic  solution  of  amplitude  decay  existed;  when  Ri-«»  it  was  a 
stable  central  point  so  there  existed  the  periodic  solution). 

The  second  order  system  was  obtained  by  expanding  the 
nonlinear  terms  and  then  retaining  the  second  order  terms.  Our 
analysis  indicated  that  the  topological  structure  was  completely 
consistent  with  the  first  order  system.  The  central  point  of  the 
first  order  system  was  still  the  central  point  of  the  second 
order  system.  The  periodic  solution  of  the  second  order  system 
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satisfied  the  famous  K^V  equation,  and  it  was  an  elliptical 
cosine  wave.  Under  other  conditions  the  solutions  of  the  second 
order  system  (which  was  the  velocity  which  varied  only  with  time) 
were  different  from  those  of  the  first  order  system.  This 
further  reflected  the  actual  situations  associated  with  many 
atmospheric  and  oceanic  phenomena  (including  turbulence). 

I .  INTRODUCTION 

The  density  of  the  atmosphere  or  the  ocean  varies  with 
height  (this  is  being  called  "layered"  or  "stratified"),  and  the 
velocity  also  varies  sharply  with  heigh't  (this  is  called  the 
"vertical  shear").  The  most  fundamental  wdve  within  the 
stratified  shear  flow  is  the  internal  gravity  wave.  The  medium 

t 

to  small  scale  atmospheric  systems  such  as  typhoon,  scowl  line, 
and  cumulus  cloud;  as  well  as  temperature  jump  and  turbulence  in 
the  ocean  are  all  closely  related  to  the  internal  gravity  wave. 

In  the  past  the  discussion. concerning  the  internal  gravity  waves 
was  mainly  concentrated  in  linear  theories.  Many  important 
conclusions  were  obtained  this  way,  such  as  the  identification  of 
the  Brunt-Vaisala  frequency  N  as  the  upper  bound  of  the  frequency 
of  internal  wave  <•>  ,  and  that  the  internal  wave  is  stable  when 
Ri  »  N^/(du/dz)^  is  greater  than  1/4^^  (here  u( z )  is  the 
fundamental  shear  flow  field).  It  was  also  established  that  Ri  - 

1/4  is  the  boundary  between  laminar  flow  and  turbulent  flow  in 

- -  2 

stably  stratified  atmosphere  (N  >  0).  In  recent  years 

observations  have  revealed  large  amplitude  fluctuations  of  the 

boundary  layer  at  night^^  and  blocked  structures'^. 

r  2] 

Intermittent  turbulent  flow  and  solitary  waves  have  also  been 
observed.  None  of  these  phenomena  can  be  explained  by  the  linear 
theory.  What  are  the  characteristics  of  nonlinear  stratified 
shear  flow?  What  kinds  of  relations  do  they  have  with  the  Ri 
number?  Can  we  find  analytical  solutions?  All  of  these  are 
topics  that  we  are  concerned  about.  Some  of  the  relatively  early 
studies  concerning  the  theoretical  examination  of  nonlinear 
Internal  waves  include  the  work  of  Benney  (1966)^^  and  Benjamin 
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(1966)^^  (1967)^^.  Recently  Maslowe  and  Redekopp*-®-^  have 

conducted  a  systematic  research  into  the  nonlinear  stratified 

shear  flows.  All  of  these  studies  applied  the  multiple  scale 

perturbation  method.  This  paper  analyzed  the  set  of  equations 

which  describes  the  flow,  expanded  the  nonlinear  terms  of  this 

[91 

set  of  equations  by  asymptotic  expansion  and  conducted  an 
analysis  on  the  phase  plane.  In  addition,  we  have  also  obtained 
the  analytical  solution  of  the  second  order  system  which  retained 
the  second  order  terms.  We  obtained  an  elliptical  cosine  wave 
where  there  was  no  shear  (du/dz  =  0).  This  paper  carefully 
examined  the  effects  of  the  shear  du/dz  and  Ri  number  on  the 
stability  of  the  flow  within  the  stratified  shear  flows.  We 
pointed  out  that  the  property  of  stability  was  different  with  a 
difference  in  the  sign  of  the  shear  du/dz.  Ri  =  1/4  was  not  a 
line  separating  stability  from  instability. 
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II.  BASIC  SET  OF  EQUATIONS 

By  using  p(z)  to  represent  density  stratification  and  u(z) 
to  represent  velocity  shear;  the  two  dimensional  (x,z)  motion  of 
a  Boussinesq  fluid  can  be  considered.  The  nonlinear  equation  of 
motion,  adiabatic  equation  and  continuity  equation  which  describe 
the  stratified  shear  flow  can  be  represented  as 


dm  .  t-  .  \  9u  ,  da  1  dp 

K) ,  r .  ,  Kf)  .*  . 

dt  dx  g 


aT+ 


v\0  *.  v'v'..: 


Here  u,w  are  the  x  and  z  components  of  velocity,  respectively,  p 
is  pressure,  p  is  density,  g  is  gravitational  acceleration, 

4v»  -  —  L  4L 
P  da 

is  a  stratification  parameter,  du/dz  is  a  velocity  shear 
parameter. 

Due  to  the  existence  of  the  nonlinear  advective  terms  in 
equation  (1), 

d(— ) 

•XT'  •¥•  — sfi 

9t  dx 

it  is  generally  difficult  to  obtain  an  analytical  solution  which 
can  explain  the  motion  phenomenon.  Here  we  will  consider  a  type 
of  progressive  wave  motion  described  below.  We  will  assume  that 
the  wave  solution  has  the  following  form: 

» -  t/(f)  *-»'($)  p-P(5) 

f  —  ti>/ 

Here  k,  n  are  wave  numbers  along  the  x,z  direction,  respectively, 

<a=kc  is  the  angular  frequency,  c  is  the  wave  speed.  None  of  the 

parameters  w,  k,  c  will  vary  with  x,z,  or  t. 

Substitute  (2)  into  (1)  we  have  _  . 

(-*  +  kU  +■  ~  w - 

6M  p 

(-»  +  W  +  &W - -  nT  -  ng 

a 

v  (3) 

(-«  +  +  *•)  ir  -  —  w  -  o 

*1/’  +  nW’  -  0  * 

Here  the  symbol  is  used  to  represent  derivative  taken  against  / 1 23 

J. 

By  eliminating  P  and  W  using  the  first  and  the  second 
equations  in  (3),  we  can  obtain  two  ordinary  differential 
equations  in  (J  and  11 . 

(«a  +  *')(—»  +  kU  +  ku)U'  —  »kgn  ■+■  nk.  ~  V 

as  (4) 

(-«  +  KU  +  W - 

tn 


*  M  *  •  *  a  *  »  *  .  *  j  *  .  •  «  *  m  *  •  *  •  *  •  >  ' 


r-  *'’*  «*•  . 


If  (-in  +  kU  +  ku)  4  0,  (4)  can  be  transformed  into 


n{gn  +  n{£Lu 
u’~— + 

-Viv 
n-7 - r^r 

(  — oj  -r  ^1/  *f  k.u) 

Here  F  (U,  tl),  G(U)  are  nonlinear  functions. 

The  right  hand  side  of  equation  (5)  does  not  contain  l> 
explicitly  so  it  is  a  set  of  ordinary  differential  equations  for 
the  self-governing  system.  Based  on  qualitative  analysis  of 
ordinary  differential  equations  it  is  7easy  to  see  that  the 

origins  u*  0,  U  =  0  is  the  singular  point  on  the  (U,n  )  phase 
plane.  At  this  point  U=0  and  p=0  so  this  point  represents  the 
undisturbed  static  condition  based  on  physical  consideration. 

In  order  to  analyze  (5),  we  will  conduct  Taylor  expansion  on 
the  nonlinear  F  and  G  terms-  on  the  right  hand  side  of  equation 
(5)  in  the  vicinity  of  the  origin  of  the  II  ,  U  plane.  We  can 
thus  obtain 


U'  — - - - V  h - 21 - JI  -  - _ — _ I/* 

(«*  +  ^)(— c  +  5)  C*1  +  k*)(— c  +  w)  (»*  +  $*)  (—  c  +  u)* 


- Zt - jju  + 

(✓  +  **)(-*+«) 


JT  —  - T"~v  U1  + 

gn(—c+*)  gn{—  e  +  m)1  t 


We  will  now  explain  that  equation  (6)  is  an  equation  which 
can  be  used  for  the  discussion  of  internal  gravity  wave.  We 
shall  place  our  emphasis  on  the  commonly  seen  practical  situation 
where  C>  u  >0.  In  order  to  simplify  our  description,  we  will 
first  discuss  the  case  where  k<  n,  and  we  will  keep  the  right 
hand  side  of  equation  (6)  up  to  the  second  order  terms.  Equation 
(6)  can  thus  be  transformed  into 


u'  —  —  u-  —  n-  —  u1-  —  nu 


- 1-(U  +  c )  u  +  en) 

nc  \  dx  / 

IT  -  —  U  +  U‘  -  ( V  +  c)V 

gne  gnc J  gnc 
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Equation  (7)  is  the  second  order  self-governing  system  which 
will  be  used  to  analyze  the  nonlinear  internal  gravity  waves  from 
now  on . 

III.  RESULTS  OF  LINEAR  WAVES 

We  will  first  discuss  the  condition  of  linear  internal 
gravity  waves  by  using  equations  (6)  and  (7).  If  we  were  to  take 
only  the  linear  terms  on  the  right  hand  side  of  equation  (6),  we 
can  use  these  terms  to  describe  the  linear  gravity  waves  at  an 
infinitesimal  distance  away  from  the  equilibrium  point  (0,0).  At 
this  time  equation  (6)  becomes 


--v' 

“  gn(—e  +  5)  U 


W  +  VK-c  +  5) 


If  there  is  only  density  stratification  p(z)  and  no  velocity 
shear  du/dz=  0  (the  same  as  Ri  -•  ),  and  we  also  assume  N  is  a 
constant,  we  can  then  obtain  the  following  from  (8) 

U "  -f - — - U  —  o  (  g ) 

If  the  variation  of  U  with  £  is  periodic,  we  can  obtain  the 
dispersion  relationship  from  (9): 


(-<•  +i)‘- 


*  +  ? 


Equation  (10)  is  exactly  the  well  known  result  of  stratified  flow 
,  and  we  can  also  come  to  the  conclusion  that  u»<N. 

If  we  now  include  the  velocity  shear  in  the  discussion 
(du/dz«0)  and  assume  that  du/dz  varies  with  z,  we  can  obtain  the 
following  equation  from  (8) 


1 

f  *  1  \ 

dm 

1 

1 

in1-)-**"  - 

da 

(-<+«) 

By  setting  k<<n,  (11)  can  be  simplified  as 


<r  u  (  du  i 

1  A-1  (  d~  - 

u "  -  -1  - ' -  y 

a1  u  —  c  («  —  £■)' 


(12) 


Here  Ri  r- 


is  the  Richardson  number. 


Equation  12  is  consistent  with  the  results  of  Taylor  and 
Goldstein^^.  Furthermore,  if  du/dz=0  equation  (12)  become  (9). 

If  we  were  to  take  only  the  linear  terms  in  equation  (7),  it 
can  be  shown  below  that  it  is  necessary  for  Ri  >1/4  in  order  to 
obtain  a  stable  wave  solution.  By  collecting  only  the  linear 
terms  in  equation  (7),  it  becomes 


v»-~u~  —  n 

nc  nc 
AT* 

IT  —  T —  U 

gnc 


(13) 


By  assuming  du/dz  and  to  be  constants,  we  -can  obtain  the 
following  equation  from  (13) 


V"  +  ii  u’  +  ^  U  -  0 

IK  »V 


We  will  set  the  U  in  equation  (14)  to  be 

it 

17  —  Ve'*"*  (15) 

We  can  then  eliminate  the  U'  in  equation  (14)  and  obtain 

{— Y 

'  -o  n A) 


4  f 


(16) 


It  is  obvious  that  it  is  necessary  for  the  coefficient  in  front 
of  V  to  be  positive  in  order  for  equation  (16)  to  have  a  wave 
solution.  This  is  to  say  that 


Ri  >  — 
4 


(17) 


All  of  the  above  are  consistent  with  the  existing  results  of 
stratified  flow  theories.  This  suggests  that  the  basic  equation 
sets  (6)  and  (7)  accurately  reflect  the  results  of  linear 
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internal  gravity  waves.  If  we  keep  the  second  order  terms  on  the 
right  hand  side  of  either  equation  (6)  or  (7),  we  can  then  use 
them  to  discuss  the  finite  amplitude  nonlinear  internal  gravity 
waves . 

IV.  THE  TOPOLOGICAL  STRUCTURE  OF  STRATIFIED  SHEAR  FLOW 
From  (13)  we  can  obtain  the  integration  curve  in  the 
vicinity  of  the  equilibrium  point  (0,0)  on  the  phase  plane  ( fl 
,U).  It  is  determined  by  the  following  equation 


du 

-nzu-gn 


The  property  of  the  phase  path  is  determined  by  the 
characteristic  root  X  of  the  characteristic  equation 

|o  —  a  £  I 


(18) 


(19) 


This  is  to  say  that 


a‘  +  — a  +AP  —  0 

dx 


(20) 


The  characteristic  root  obtained  from  (20)  is 


du  du  , - 

—27  ± -27V 


(21) 


This  is  why  the  phase  paths  (n  ,U)  in  the  vicinity  of  the 
origin  (0,0)  of  the  parameter  (du/dz,N  )  are  divided  into  several 
different  areas  as  shown  in  Figure  1.  The  black  dot  •  in  this 
figure  represents  the  equilibrium  point  which  is  stable,  the 
white  dot  o  represents  the  unstable  equilibrium  point,  and  the 
direction  of  the  row  indicates  the  direction  in  which  time  is 
increasing.  Since  £  =  kx  +  nz  -  wt,  so  for  a  certain  given 
point  (X,Z)  the  direction  in  which  time  is  increasing  is  the 
direction  in  which  £  is  decreasing.  The  equation  for  the 


dashed  line  in  this  figure  is 


(V)“T  (t^’  (c) 

4  dz  ■ 

It  is  a  parabola. 


Figure  1 

1.  stable  nodal  point  area 

2.  stable  focal  point  area 

3.  unstable  focal  point  area 

4.  unstable  nodal  point  area 

5.  saddle  point  area 

6.  unstable 

_  2 

It  can  be  seen  from  this  figure  that  the  parameter  plane  (du/dz, N) 
divides  the  properties  of  the  solution  of  the  stratified  shear 
flow  (the  phase  diagram  which  used  11  as  abscissa  and  U  as 
ordinate)  into  five  areas  (i),  (ii),  (iii),  (iv),  (v)  and  four 
boundaries  (vi),  (vii),  (viii),  and  (ix): 

(i)  Ri<  0  unstable  saddle  point  area;  x,  >0,  xa  <0 

(ii)  Ri<  1/4,  du/dz  <  0  stable  nodal  point  area;  xa  >xa  >0 

(iii)  Ri>  1/4,  du/dz  <  0  stable  focal  point  area;  xa  and  xa 

are  complex  conjugates  with  real  parts  being  positive. 

(iv)  Ri>  1/4,  du/dz  >  0  unstable  focal  point  area;  and  Xa 
are  complex  conjugates  with  real  parts  being  negative. 


(v)  Ri  <1/4,  du/dz  >  0  unstable  nodal  point  area;  x2 <  X, <0 

(vi)  Ri-«,  du/dz=  0  stable  central  point;  x2  and  x2  are 
pure  imaginary  roots. 

(vii)  Ri  =  1/4,  du/dz  <0;  xl  =  X2  >  0 
(viii)  Ri  =  1/4,  du/dz  >0;  X2  =  X2  <  0 

(ix)  Ri  =  0;  X2  is  either  positive  or  negative, x2  =0 
It  should  be  pointed  out  that  the  above  analysis  was  done 
under  the  conditions  of  C>  u,  k<<  n. 

For  the  condition  where  c<  u  the  above  analysis  can  still  be 
applied. 

All  we  have  to  do  is  to  switch  the  du/dz>  0  (or  du/dz<  0)  with 

du/dz  <0  (or  du/dz  >0).  This  is  to  say  that  the  first  quadrant 

_  2 

of  the  parameter  plane  (du/dz, N  )  is  a  stable  area  while  the 

second  quadrant  is  an  unstable  area. 

For  the  condition  where  k  is  not  smaller  than  n,  all  we  have 

n*  +  k2 

to  do  is  to  replace  the  parameter  Ri  with  - —  Ri  in  the  above 

analysis.  This  is  to  say  that  the  stable  area  is  still  related 

-  • 

to  k  and  n. 

In  this  way,  we  can  classify  the  flows  according  to  their  Ri 
numbers. 

Up  to  now  we  have  only  analyzed  the  first  order  system  (13), 
now  we  will  analyze  the  second  order  system  (  7).  According  to 
the  theory  of  Poincare-Bendixson^  ,  since  the  nonlinear  terms 
in  (7)  are  separately 

X(il,  w 

gncf 

du_  (22) 

y(n,to  —  lLu,--*-nu 

tie*  me 1 

so  we  will  have 

Wn,  ft)  -  n)  -  gjL(0>_0)  _  9Y(0,  0)  _  9X(0,  0)  _  07(0,  0)  _  p 

du  au  an  an 

In  this  way  the  qualitative  properties  of  the  second  order  system 
(  7)  in  the  vicinity  of  the  origin  (0,0)  under  the  conditions  of 
(i),  (ii),  (iii),  (iv),  (v),  (vi),  (vii),  and  (viii)  are 
completely  the  same  as  the  analysis  done  on  the  first  order 


system  (13). 

At  the  same  time,  relative  to  the  central  point  (Ri-«,  du/dz 
=  0)  of  the  first  order  system,  the  second  order  system  (7) 
becomes 


tr _ —  n  -S-nv  -  Bin,  u) 

nc  nc 

nr  -  —  v  +  -^T  V2  -  Ain,  u) 

gnc  gf>c 


(24) 


[13] 

According  to  the  symmetry  principle  ,  we  have 


.an,u)-  A{-n,u)  Bin,  £/)-  -  Bi-n,  u)  (25) 

We  can  be  sure  that  the  central  point  of  the  first  order  system 
is  still  the  central  point  of  the  second  order  system.  Even 
under  condition  (vi),  the  qualitative  property  of  the  second 
order  system  is  the  same  as  that  of  the  first  order  system. 

This  is  why  the  qualitative  structure  of  the  solution  of  the 
first  order  system  (13)  is  completely  the  same  as  that  of  the 
second  order  system  (7)  in  the  vicinity  of  the  origin.  The  only 
difference  is  that  the  rate  of  variation  of  U  and  fl  with  t  (or 
5) is  different  on  the  phase  diagram. 

From  the  complete  phase  diagram  of  the  stratified  shear  flow 

we  can  see  that  when  C  >u ,  we  will  have: 

2 

(1)  N  <  0  (the  same  as  dp/dz>  0)  is  always  a  factor  for 
instability.  As  for  whether  the  shear  is  a  factor  for  stability 
or  a  factor  for  instability  depends  on  whether  du/dz<  0  or  du/dz> 
0.  In  the  past,  the  discussion  of  a  Kevin-Helmholz  wave  always 
treated  the  shear  as  a  factor  for  instability. 

(2)  Miles  (1961)  and  Howard  (1961)  gave  the  necessary 
condition  for  a  stable  internal  gravity  wave  as  Ri>  1/4,  and  they 
did  not  include  the  condition  where  0<  Ri<  1/4.  It  can  be  seen 
from  the  phase  diagram  that  as  long  as  du/dz<  0,  the  wave  is 
still  stable  when  0  <Ri<1/4. 

(3)  In  the  past,  people  always  believed  that  the  wave  is 
always  stable  when  the  Ri  number  is  very  large  (at  this  time  the 
shear  du/dz  is  very  small).  This  is  why  the  production  of 

o 

turbulence  under  the  condition  of  stable  stratification  (N  >  0,  dp 
/dz  <0)  could  not  be  explained.  According  to  our  analysis, 
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instability  can  still  be  reached  under  very  high  Ri  as  long  as 
du/dz  >0. 

(4)  When  1/4  <Ri<®  ,  there  exist  periodic  solutions  with 
either  amplitude  growth  (under  the  condition  of  du/dz>  0)  or 
amplitude  decay  (under  the  condition  of  du/dz  <0).  When  Ri-®  , 
there  is  a  periodic  solution  and  when  Ri  ^1/4  there  is  no 
periodic  solution.  During  the  night  there  is  a  phenomenon  of 
frequent  appearance  or  vanishing  of  large  amplitude  waves  in  the 
boundary  layer^].  We  believe  that  this  is  exactly  the  condition 
of  the  stratified  shear  flow  within  the  focal  point  area  where 
1/4  <Ri  <-.  r 

V.  THE  NONLINEAR  SOLUTION  OF  THE  STRATIFIED  SHEAR  FLOW 

In  order  to  compare  the  rates  of  variation  of  U  and  *  with  $ 
(or  t)  between  the  nonlinear  system  (second  order  system)  and  the 
linear  system  (first  order  system),  we  will  now  derive  the 
analytical  solutions  U  (  tj-  and  11(0  for  the  second  order  system 
(7). 

By  taking  the  derivative  against  5  ,  the  first  equation  in 
(7)  becomes 

^  ^  26 ) 

ir  -  —  —  tr  -  -2-  it  —  ♦  zau’  -  —An’u  +  irn) 

m  w  *c  *r 

We  can  now  obtain  the  second  order  nonlinear  equation  of  U  by 
substituting  the  in  the  second  equation  of  (7)  and  the  n  '  in 
the  first  equation  (7)  into  the  formulation. 

°"  +  5(1  +  t)'''-T77^  +  |:1'(1+7)'-0  <”> 

C 

We  will  now  derive  the  analytical  solution  of  the  nonlinear 
equation  (27)  as  follows.  To  serve  as  a  contrast,  we  will  first 
consider  a  relatively  simple  condition.  This  is  to  say  that  we 
will  first  eliminate  the  nonlinear  terms  from  equation  (27)  and 
then  try  to  solve  the  linearized  equation: 


V”  +  —  U'  U  -  0 
ne  nle 
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(28) 


Equation  (28)  is  actually  equation  (14),  so  its 
corresponding  characteristic  equation  is 


r*  +  —  r  +  If- 1  -  0 
nc  »f 


The  characteristic  root  of  (29)  is 


i5  dir 
da  ,  da 


——  ± 


(29) 


(30) 


Here  x  is  the  representative  equation  (21). 

The  form  of  the  solution  of  equation  (28)  U  (5)  can  thus  be 
determined  based  on  whether  the  characteristic  root  r  (or  x )  is  a 
real  root,  an  imaginary  root  or  a  complex  conjugate  root.  We  can 
then  obtain  11  (5  )  from  the  second  equation  of  (13).  These  special 
solutions  are  listed  as  follows: 

(1)  When  x  is  a  real  root,  the  U  (5 )  and  * (5 )  corresponding 
to  the  two  straight  lines  which  pass  through  the  saddle  point  on 
ncdal  point  along  the  phase  path  are 


17(f)  -  e-t*  _ 

11(f)  -  e-t*  _  ? JL 

fi 


(31) 


(2)  When  x  are  complex  conjugate  roots,  the  solutions  of  U(5) 
and  * (C )  corresponding  to  the  spiral  line  in  the  focal  point 
area  on  the  phase  diagram  are: 


J  -  «**  *eoa^-  {-  «••  (^x  +  na  —  o«> 

w 

n<£i- it  •  «»(£  f  +  - - *  my  * 


(32) 


Here  <*  and  &  are  the  real  and  Imaginary  part  of  x,  respectively; 
and  tgA  »  -  VJrTTT 

(3)  When  X  is  a  pure  imaginary  root,  the  U(£)  and  11(0 
corresponding  to  the  .closed  trajectory  around  the  central  point 
of  the  phase  diagram  are 
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'•  A  .**  /■ 


*.  \  / 


•  V 


/  J 


l/({)  ™  eo*  —  f  —  co«  —  (^x  +  «*  —  u/) 

»<•  BC 

n(f )  “  “  7  “®  “  f  “  —  y  * in  —  (^*  +  **  —  mi) 


(33) 


We  will  now  solve  for  the  nonlinear  solutions  corresponding 
to  equations  (31)  and  (33)  as  follows.  Since  the  relationship 
between  U  and  fl  is  given  by  the  phase  diagram,  we  will  not  try  to 
fintf  the  solution  by  directly  using  equation  (27).  Instead,  we 
will  use  the  original  set  of  equations  (7). 

For  the  two  straight  lines  passing  through  the  saddle  point 
or  the  nodal  point  on  the  phase  diagram,  we  can  derive  the 
following  equation  from  (31) 


V  —  —  n 


By  substituting  (34)  into  (7)  we  can  obtain 


By  integrating  (35)  we  can  get 


£L  -  i  i/(t/  +  c) 


(34) 


(35) 


v 

U  +  c 


(36) 


(37) 


£  » 


So  we  set 


U  -- f(l  +  --  j-  [l  +  +  »»-  «*)] 


(38) 


By  substituting  (38)  into  the  second  equation  of  (7)  we  can  get 


n"  (‘ I1 +  »*-«'>] 


(39) 


'.s  v. 


15 


By  comparing  (38)  and  (39)  with  the  linear  solution  (31),  we  can 
see  that  the  velocities  at  which  they  approach  (or  get  away  from) 
the  origin  (the  equilibrium  point)  are  different.  The  solution 
obtained  by  considering  only  the  linear  terms  has  the  form  of  an 
exponential  function.  The  solution  obtained  by  considering  the 
nonlinear  terms  has  the  form  of  a  hyperbolic  tangent  function. 

What's  more  meaningful  is  that  (38)  and  (39)  are  very 
similar  to  the  progressive  wave  solution  of  the  nonlinear 
dispersion  equation  (which  is  the  Burgers  equation) 


_  a** 
a,  +“a 7  “a? 


(v  is  the  dispersion  coefficient) (40) 


The  progressive  wave  solution  with  u  >  U(g)  and  5  *  kx-w  t  is 


V 


t 


(41) 


In  equations  (38)  and  (39) 

1 

-jjp  is  equivalent  to  u  (42) 

So  the  IM  value  increases  as  du/dz  increases,  corresponding  to 
decreasing  u  value.  In  this  way,  U  will  vary  sharply  at  the 
location  where  5  changes  sign.  This  explains  that  for  the 
nonlinear  solutions  of  (38)  and  (39)  there  is  a  phenomenon  of 
sudden  variations  in  U  and  U  when  du/dz  is  very  large.  This 
phenomenon  will  not  occur  if  we  only  consider  the  linear  system. 
The  conclusion  of  this  paper  can  be  used  to  explain  a  lot  of  the 
intermittent  phenomena  in  the  atmosphere  on  the  ocean. 

Now  we  will  derive  the  nonlinear  solution  corresponding  to 
(33)>  at  the  time  the  origin  of  the  phase  plane  is  the 
central  point.  From  (18)  or  (33)  or  (24)  we  can  obtain 

dU _ jn_ 

*u  ¥u  (43) 

f 

.1 

So  the  closed  phase  path  around  the  central  point  is 


n *  +  2!  -  d 


(  Ah  \ 


In  this  equation  D  is  an  integration  constant.  (44)  is  an 
ellipse. 

By  taking  the  derivative  against  5  for  the  first  equation  of 
(24), then  substituting cne  second  equation  of  (24)  and  equation 
(44)  into  the  formulation,  and  retaining  only  up  to  the  second  order 
terms,  we  can  get 


Un  —  —  —  i"  —  —  o  +  1 


(45) 


By  taking  the  derivative  of  (45)  with  respect  to  £  we  can  obtain 
the  famous  K^V  equation  7 


IT*  +  W  +  ~  V’  -  0 


(46) 


We  have  already  obtained  the  solution  of  (45)  and  (46),  it  is  an 
elliptical  cosine  function^^cn 

u(s )  -  Vi  +  (i/,  -  ~  u*>  f  **  "  w,)  (47) 

Here  >0,  U2  <0,  U2<U2<0  is  a  third  order  equation 

Him-W  +  +  B  (48) 

These  are  independent  single  real  roots,  and  B  is  a  constant. 

It  can  be  seen  from  the  comparison  between  (47)  and  the 
corresponding  linear  solution  (33)  that  the  linear  wave  has  the 
form  of  a  cosine  (or  sine)  function  while  the  solution  obtained 
by  taking  into  consideration  the  nonlinear  effects  is  an 
elliptical  cosine  function.  The  periodicity  and  wavelength  are 
also  different.  Furthermore,  the  amplitude  of  a  nonlinear  wave  is 
related  to  its  wave  speed*'  ^ . 

Now  we  have  obtained  nonlinear  solutions  for  the  two 
straight  lines  passing  through  the  saddle  point  or  the  nodal 
point  of  the  phase  diagram  as  well  as  the  closed  curve  around  the 
central  point  of  the  phase  diagram  (  f,U). 


VI.  CONCLUSION 

We  have  obtained  the  structure  of  geometric  topology  for  the 

phase  diagram  of  the  stratified  shear  flow  in  the  vicinity  of  the 

origin  of  the  (n,U)  phase  plane.  Its  property  is  determined  by 

2 

the  stratification  parameter  N  and  the  shear  parameter  du/dz. 

Our  analysis  indicates  that  the  stability  of  the  flow  depends  not 

2 

only  on  the  sign  of  stratification  N  (same  as  3p/9z)  but  also 

on  the  sign  of  the  shear  (du/dz).  This  is  why  only  the  first 

quadrant  (c<u)  and  second  quadrant  (c>u)  on  the  parameter  plane 
_  O 

(du/dz,  N  )  are  stable.  Ri  «  1/4  is  not  the  boundary  between 
stable  and  unstable.  It  is  actually  the  boundary  between  periodic 
and  nonperiodic  solutions. 

The  qualitative  properties  of  a  linear  system  are  completely 
the  same  as  a  nonlinear  system,  but  the  rates  of  variations  are 
different.  For  the  nonlinear  system  there  is  a  phenomenon  of 
rapid  variation  when  ldu/dz.1*  is  very  large.  The  amplitude  of  a 
nonlinear  periodic  wave  is  related  to  its  wave. speed. 

The  authors  have  received  much  guidance  and  assistance  from 
Professor  Huang  Dun  and  Associate  Professor  Liao  Ke  Ren  of  the 
Dpeartment  of  Mathematics,  Beijing  University.  We  would  like  to 
express  our  sincere  appreciation. 
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